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ABSTRACT 

This study reviews the theories and modelling methods for describing interfacial 
delamination failure process between two bonded cementitious materials. Complex interfacial stress 
conditions at discontinuities and areas of high stress concentrations were primary areas of concern. 
Distinct analytical cases involving intrinsic material and structural property variables were 
considered. An approach based on plane strain analysis within the context of Interface Cohesive 
Zone Model (ICZM) was cited and presented as viable for describing and predicting delamination 
mode of failure in bonded concrete overlays systems (BCOs). The study shows that the use of 
numerical computational tools is vital in resolving the manifold complexity associated with 
interfacial delamination problems. In the concluding analytical model, it is evident that the numerical 
values of the delamination failure coefficient(D) and the corresponding Mixed-Mode energy release 
rates (G ic ) vary depending on the overlay structural scale, the type of problem (plane stress or plane 
strain) and the degree of mismatched properties between the overlay and the substrate. 

Keyword: Interfacial, ICZM, Delamination, BCOs, Mismatched. 

1.0 INTRODUCTION 

Adequate interfacial bond performance of Bonded Concrete Overlays (BCOs) requires novel 
integration of material mixture design, compatibility model development, and robust interfacial 
bonding techniques. This whole process entails the use of the right material, on the right substrate, in 
the right way, in order to secure the best possible composite behaviour. In this respect, the structural 
integrity of pavement can be reinstated with enormous benefits, ranging from resource conservation 
to good returns on investment. 

However, in spite of the plausible benefits accruing from BCO system of repair when 
compared to Un-bonded Concrete Overlay systems (UBCOs), early-age delamination problem 
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remains a bugging issue (Karadelis and Koutselas, 2003; Olubanwo, 2013). Hence, the doubt often 
arises as to whether the bond integrity of the BCO systems can be trusted, in particular at early-age, 
when the interface experiences strident self-equilibrated stresses resulting from both mismatched 
properties and differential length change effects. 

Certainly, solving delamination problem requires a more pragmatic approach than just 
ensuring early interfacial bond development. Optimum solution often requires a holistic approach 
(Morgan, 1996; Emmons and Vaysburd, 1996) during which the thermo-mechanical compatibility 
and stability between the bonded layers are assessed and ascertained a priori. Implementing such 
tasks is non-trivial. Typically, it may involve bond optimization analysis using empirical and 
computational methods in order to assess and predict both early and long-term durability 
performance of the bonded system. Extensive review and work in this line has been presented 
elsewhere (Olubanwo and Karadelis, 2014). Until now, much of the optimization works for BCO 
system design revolves largely round thickness requirement. Fragmentary approach of this nature is 
somewhat insufficient and could result in abrupt material and subsequent structural failure of the 
BCOs (Olubanwo, 2013). In the literature, several successful interface models exist, not without 
drawbacks though. The following sections review and discuss briefly some prominent ones and their 
governing theories. 

2.0 DELAMINATION MODEL DESCRIPTIONS 

From experimental and analytical standpoints, two basic approaches are commonly used for 
simulating and describing the failure or de-bonding process of the interface- (1) stress-basedfailure 
criterion approach and (2) Energy-based fracture criterion approach. By definition, the two 
approaches defer in experimental concepts and computational techniques. For instance, in the limit 
analysis or the so-called stress-based approach, the interface is often assumed as perfectly bonded, 
while the classical energy-based method (Linear Elastic Fracture Mechanics - LEFM)treats the 
interface as having some well-established intrinsic defects. 

However in quasi-brittle materials like concrete, the two extreme collapse methods cited 
above are rear and mostly infeasible; hence, more robust methods are desirable. In this paper, 
Interface Cohesive Zone Model (ICZM) within the concept of non-linear fracture mechanics is 
reviewed and proposed as a desirable alternative for describing the failure mechanism of concrete 
based interface. In the method, delamination process involving both crack initiation and propagation 
within a unified model is represented. The model treats delamination as progressive. Special cases of 
failure due to differential length change at discontinuities involving possible elastic mismatched 
conditions and structural scale properties are considered distinctly. 

2.1 Traditional Stress-based approach 

Essentially in stress-based approach, in order to ensure that the delamination process of the 
interface is adequately depicted, the mechanical characterization of the interface requires two basic 
descriptive states: (1) a state representing a perfectly bonded condition, and (2) a state defining 
delamination on set and propagation. In the former, adhesion between the bonded layers is assumed 
sufficiently strong; there by, imposing both stress and displacement continuity across the interface. 
In Shah and Stang (1996), the corresponding kinetics and the kinematics for continuity requirements 
at the interface are given by: 




(2.1) 



86 



International Journal of Civil Engineering and Technology (IJCIET), ISSN 0976 - 6308 (Print), 
ISSN 0976 - 6316(Online), Volume 6, Issue 5, May (2015), pp. 85-99 © IAEME 



WhereJbisall points on the bonded interface, x and c are respective shear and normal stresses 
of a point on the interface, while v and u are their corresponding tangential and normal 
displacements. Note, the subscripts 'top' and 'btm' represent the top and bottom layer respectively. 

As seen, equation 2. 1 lacks practicality in many instances, basically because it is premised on 
the assumption of a perfectly elastically bonded interface with no possibility of yielding or de- 
bonding. In reality, the interface yields or de-bonds at considerable lower stresses compared to 
adjacent bulk domains as stresses localize or concentrate in the plane of the interface due mainly to 
bond imperfection and mismatched elastic properties between the bonded materials. It has been 
shown that such localized or concentrated stresses can be three times more detrimental than the 
average stresses developing in the adjacent bonded materials (Kirsch, 1898; Dantu, 1958). Hence, 
characterizing the interface with a finite strength is commonplace in practice. Typically, for 
composite interface model of this nature, a general stress-based failure criterion takes the form: 

F(T,a,Pj) = 0 i=l....,n (2.2) 

Where, Pjis one of n strength parameters, while other parameters are as given in equation 2.1. 

From equation 2.2, the interfacial de-bonding process in limit analysis now permits the 
bonded interface to separate once it is loaded beyond its critical bond strength. The governing 
constitutive relations in this case are generally based on the kinematics of the interface, as the 
interface changes from its continuity condition to a prescribed surface traction boundary condition. 
In this case, the failure criterion takes the form (Shah and Stang, 1996): 



x top x btm f 
v top = v btm — g 
U t0 p - U btm > 0 



• on I d (2.3) 



Where, I d is all points on the de-bonded interface, while T and 'g' are prescribed surface 
tractions in the general case. All other parameters are as given in equation 2.1. 

Though the conditions given in equation (2.3) show a high explicit level about the nature of 
the de -bonding, the possibility of surface overlapping during de-bonding process is precluded; and 
often, this can be difficult to substantiate in reality, especially in cementitious materials where 
complicated interfacial contact problems dominate during de-bonding initiation stage (Shah and 
Stang, 1996). During subsequent steps, and at critical cracking stage of the interface, the interface 
attains a stress-free state (Atkinsonet al., 1982; Stang and Shah, 1986; and Morrison et al., 1988); 
hence, the frictional stress givenin equation (2.3) vanishes accordingly, so that it reads: 

x top = x btm = 0) T * \ 

v _ v _n onI d (2.4) 

v top - v btm - U J 

If the interface is simultaneously influenced by stress and displacement continuities 
perpendicular to the interface, the boundary conditions given in equation (2.4) can be extended to 
include: 

°top = a btm = 0} 

„ _„ _ n onI d (2.5) 

u top — u btm — U J 

Thus far, from engineering standpoint, the use of limit analysis as demonstrated above for 
both perfectly bonded and de-bonded interface surface characterization seems reasonable and 
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acceptable, but its inability to explain or capture the infinite stress condition at the crack-tip between 
the perfectly bonded region and the de-bonded zone of the interface remains a major drawback. The 
analytical proof of such stress singularities are based on complete linear elastic solution of the de- 
bonded interface problem. In the literature, it has been shown that the presence of these singularities 
is the reason why finite element analysis involving interface problem is mesh-dependent 
(Mormonieret al., 1988). In contrast, the justification of infinite stress state is also hard, if not 
impossible, if it is generally accepted that no material can be loaded beyond its yield strength. 
Consequently, it is inferred that many of the parameters employed in determining stress-based de- 
bonding criterion cannot be regarded as absolute material parameters, knowing that their values vary 
widely according to the type and complexity of the test and analysis used. Besides, it is likely - due 
to the incompleteness of the analysis - that the value of the strength parameters corresponding to a 
given interface will depend on the size, geometry and loading conditions of the composite system. 
Size effects, for instance, have been observed in similar tests, which, for instance, pure shear analysis 
fails to predict (Shah and Stang, 1996; Bazant and Zi, 2003). With the advent of modern 
computational tools, however, many of the problems associated with incomplete analysis or rigorous 
analytical solutions for de-bonded interface can be resolved. 

2.2 Energy-based criterion 

On the other hand, in energy-based criterion, methods based on classical Linear Elastic 
Interface Fracture Mechanics (LEEFM) have been found effective in describing and modelling crack 
propagation, particularly where material nonlinearities are negligible (Turon, et. al, 2004). As 
illustrated in Figure 2.1, for a bonded bi-material interface whose adjacent domains are characterized 
linearly by Ej, u^and vj for Young's Modulus, shear modulus, and Poisson's ratio of each domain 
respectively, it has been shown that there exists an intrinsic singularity with oscillatory field ahead of 
the crack tip (Williams, 1959), caused by the asymmetry in the elastic properties across the interface. 
This oscillatory field controls the measure of the competing or complex stress state near the crack- 
tip, which in terms of stress intensity factors can be expressed as: 




(2.6) 



Material #2 

E7 , v?. 



Figure 2.1: Linear crack along a Bi-material Interface 

Wherer 18 = exp(is log r) = cos (slogr) + i sin (slogr),i = V— 1, r is the distance ahead of the 
crack-tip, s is the oscillatory index defined later in equation (2.9), K x and K 2 are components of the 
complex stress intensity factor derived by Rice and Sih (1965) by solving the following logarithmic 
expressions which result from a full boundary-value problem of a given test specimen (Carlsson and 
Prasad, 1993; Chandra, 2002): 
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v _ o[cos(slog 2a)+2s sin (8log2a)]+{i[sin (£log2a)-2£ cos (Elog2a)]} r— „. 

K x — Va {2. 1) 

COSllTO 

_ x[cos(£log2a)+2ssin (elog2a)]-{a[sin (Elog2a)-2ecos (Elog2a)]} i— 
K 2 — 7 Va (^-o) 

cosh TO 

From where, sisestimated as: 

6 = -Inp=£| (2-9) 

In equation (2.9), (/?)relates to one of Dundur's elastic mismatched parameters (Dundur, 
1969) which measures the relative compressibility of the two bonded materials, commonly estimated 
from equation (2.10), say, for plane strain problems (Mei et. al, 2007); while its counterpart (a) given 
in equation (2.11) measures the corresponding relative stiffness (Mei et. al, 2007; Schmauder, 
1990;Bower, 2010). 



2 



Mi(l-2v 2 )-M 2 (l-2Vi) 



v 2 )+M2(l + Vi) 

Which on simplifying yields: 

„ _ £ I 1 (l-v 1 )(l-2v 2 )-ir 2 (l-v 2 )(l-2i; 1 ) 



(2.10b) 



r 2(l-V 1 )(l-v 2 )(£ 1 + £ 2 ) 

Where,^- = Ej(l — vf) plain strain Young's Modulus for material i 
a=%^ (2.11) 

E ± + E 2 

Subsequently, under a Mixed-Mode fracture analysis, the energy release rate,(G), for crack 
extension per unit length along the interface for plain strain is generally given by (Carlsson and 
Prasad, 1993): 

G = — (2 12) 

E*cosh 2 nE y ' ' 

Where, 

\K\ = V^ 2 + Kl (2.13) 
cosh 2 ns = 1/(1- /? 2 ) (2.14) 

^ + (2-15) 

E* 2 E 2 J v ' 

Thus, by Mode-Mixity, the value of(G) as a function of the loading phase angle Q¥~) follows 
the real and imaginary stress intensity factors of the remote field lying ahead of the crack tip. This 
phase angle is typically expressed as: 



V = tan" 1 
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Where, L is the arbitrary reference length selected to characterize the remote field. For most bi- 
material systems, it is clear that the value and the effect of nonzero (/?) is small, and so in significant 
(Bower, 2010; Buyukozturk and Hearing 1998). Thus, by settings = 0, for most material 
combinations, equation 2.16asimplifies to: 



The corresponding displacement components behind the crack tip are given in Bower (2010) by: 



From the above equations, it is clear that the asymptotic solution for the interface crack 
differs significantly from the corresponding solution for a homogenous solid, because the oscillatory 
character due to both stresses and displacements increases frequency as crack tip is approached; 
hence, making it difficult to discretize the remote loading. Besides, the crack planes are predicted 
overlapping near the crack-tip a priori, which perhaps is still less than clear in many practical 
instances (Bower, 2010). 

As seen above, the application of LEIFM approach is attractive when considering crack 
propagation process particularly for brittle materials. Here, the critical fracture condition is assumed 
to have been reached when the energy release rate G equals the fracture toughness of the interface 
G(¥); that is: 



In many experimental instances, it has been shown (Suo and Hutchinson, 1989; 
Charalambideset. al., 1990) that this interface resistance to delamination increases rapidly with phase 
angle. 

In essence, while both approaches described above provide some degree of analytical 
comfort; in reality, evidence of initial perfect interfacial bonding or the presence of initial interfacial 
crack, together with its location and size may be difficult to spot or substantiate in a cementitious 
bonded overlay composite system. It is therefore thinkable to seek an enhanced method, where both 
interfacial crack initiation and propagation processes are described within a unified model. 
Employing nonlinear Interface Cohesive Zone Models (ICZM) affords a common opportunity for 
simulating both interface crack nucleation and crack growth. The governing concepts are well- 
known and particularly suitable for representing adhesion and decohesion processes between 
dissimilar materials (Mei et al, 2010). 

2.3 Interface Cohesive Zone Model (ICZM) 

In the Interface Cohesive Zone Model (ICZM), the primary consequence of nonlinear 
fracture analysis is based on the assumption of a finite fracture (cohesive) zone existing in the 
vicinity and ahead of the crack-tip, following Dugdale (1960) and Barenblatt (1962) models. The 
models as depicted in Figure 2.2 show that the so-called stress singularity (infinite stress state) 
concept commonly associated with crack-tips in elasticity theory is unrealistic (Cornec, et al., 2003). 




(2.16b) 




(2.17) 



G = GQP) 



(2.18) 
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Figure 2.2: Comparison between (a) Dugdale and Barenblatt Models and (b) Stress singularity in 

Elasticity theory 

The cohesive interface models represent a condition of a perfectly bonded interface 
undergoing de-bonding when the cohesive strength of the bonded interface varnishes with 
displacement discontinuity (Shah and Stang, 1996). Hence, for phenomenological nucleation and 
propagation descriptions, ICZM is most attractive; in particular in materials like concrete, where the 
non-linearity and the damage zone in the vicinity and ahead of the crack-tip cannot be neglected. 

Typically, for a complete interface nonlinear fracture model description, the following are essential: 

• The behaviour of the bulk materials, and 

• The behaviour of the fracture zone, where conditions for crack formation and evolution are 
pre-defined along a known crack path. 

In general, cementitious materials are classified as quasi-brittle. Besides, because cementitious 
materials exhibit little or no bulk dissipation, an isotropic linear elastic behaviour, characterized by 
elastic modulus (E) and Poisson's ratio (v), can be assumed. With respect to the fracture zone 
behaviour, a softening damage characteristic along the interface can be assumed based on the 
kinetics and kinematics relations defined in Figure 2.3c&d. 
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Figure 2.3: Interface configurations with Fracture Process Zone (FPZ) and Interface Stress 

Distribution 
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As seen, Figure 2.3 shows a BCO in its un-deformed and deformed configurations with a 
visible or true crack. The resulting FPZ, interface stress / conjugate variables, stress distribution 
curve between the FPZ and the elastic bi-material interface, and the constitutive relation defining 
each zone in the cohesive model after deformation are equally depicted. The overall cohesive model 
illustrated here follows the assumptions given below: 

• The FPZ along the interface localizes into a single line ahead of the crack tip, with no 
possibility of kinking. 

• The FPZ or the fictitious crack length is assumed mostly dominated by inelastic deformation. 

• The materials lying adjacent the fictitious crack behave linearly elastic. 

• The constitutive law governing the inelastic deformation at the FPZ assumes stress- 
displacement relationship. 

In this respect, the associated kinematics in this sense refers to the relative motion of the two 
deformed layers at the interface and may be described as: 

• In Figure 2.3(c), u represents a unit vector in normal direction to the interface. 

• At the interface, two mutual tangential unit vectors,vand w,are introduced. 

• But for cementitious inter face where isotropic condition applies, the tangential deformation 
along v and w directions is treated as equal. Hence, the constitutive relation can be defined in 
terms of scalar Cartesian components Au and At only (Bower, 2010), which represent the 
relative displacements of two initially coincident points at the interface, in normal and 

tangential directions respectively. (Where: At = Vv 2+ w 2 ). 

The kinetics relate to the forces acting between the two contacting layers. In this case, the two 
equal and opposite tractions are assumed acting on two initially coinciding points before and during 
interface deformation. Thus, under isotropic condition, the corresponding interface tractions are 
given by the scalar components T n and T t in thenormal and tangential plane respectively. 

3.0 APPLICATIONS OF INTERFACE COHESIVE ZONE MODEL TO CONCRETE 

The use of finite cohesive zone for cementitious materials is well-known, following the linear 
softening model of Hillerborget al. (1976). The application of cohesive zone model (CZM) for 
cementitious materials has since grown into popularity due to its computational convenience, and it 
is probably the best fracture model for simulating fracture processes in cementitious materials and 
structures (Bazant, et al., 2002). 

Interestingly, for concrete, both linear softening model, introduced by Hillerborget al. (1976), 
and bilinear softening model, developed by Petersson (1981), can be implemented. With most Finite 
Elementcodes, the surface traction can easily be obtained as an extrapolation of standard Gauss nodal 
stresses between adjacent continuum interface elements. 

For instance, in ANSYS FE software, the contacting interface can be modelled as a zero- 
thickness contact plane characterized with associative cohesive elements with constitutive properties 
along a pre-defined interface between two adjacent continuum elements and the resulting interfacial 
damage initiation and evolution defined by the nonlinear traction-separation (bi-linear) law described 
in Alfano and Crisfield(2001). 

As shown in Figure 3.1, the interface in a bi-linear model is assumed to behave elastically 
under deformation with initial stiffness (k 0 ) until the applied stress reaches the cohesive strength (<r 0 ) 
of the interface, at which point the damage initiation occurs. It should be noted that the initial elastic 
stiffness (k 0 ) has a character of a penalty factor only rather than a physical stiffness, hence it is 
discretionary kept high to ensure minimum elastic deformations of the interface (i.e. of negligible 
degree), so as to minimise interpenetration, separation or sliding prior to cracking. 
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<T(u) 




crack opening (u) 



(a) (b) 

Figure 3.1: a) Definition of stress and conjugate variables, and b) Bilinear softening relation 



As seen in Figure 3.1, the delamination process is defined by two slopes OA and AC. The 
initial slope OA represents the linear elastic regime of the curve, while the second slope defines the 
softening part of the curve in a linear function. De-bonding is assumed to initiate at peak contact 
stress (er 0 ) at point A, and grows linearly as a function of de-bonding parameter (d). The value of (d) 
evolves progressively from 0 to lbased on the conditions shown in equation 3 . 1 till all the interface 
stresses reduce to zero at critical crack point C (it£ ). 



0 for u n = u n 

0 < d < 1 for u n > u n 



Where, 



u n = separation of the interface elements over the entire loading history. 
— ff o 

u n = — = critical separation for damage initiation 
k 0 

er 0 = cohesive strength 

k 0 = initial elastic contact stiffness 



Thus, for each mode of failure during loading, the fracture cohesive stress (T) can be related 
to the opening or sliding displacement linearly by: 

T _ (a = k n u n (l - d n ) for Mode I 

{ t= k t u t (l - d t )for Mode II 

Where, a and t are the cohesive stresses in the normal and tangential directions respectively, 
while/c n and k t denote the corresponding contact stiffnesses. it n andu t represent the accompany 
displacements after deformation,while d n and d t are the resulting de-bonding parameters in Mode I 
and Mode II respectively. 

However, for bonded dissimilar materials, Mixed-Mode delamination is common during 
loading and failure process, thus, the criteria for damage initiation and final failure must account for 
the concomitant effects of Mode I and Mode II. In that respect, both normal and tangential traction- 
separation curves can be expanded in the (u n and u t ) - plane, as illustrated in Figure 3.2.From here, 
it is clear that the normal and shear stresses depend not only on their corresponding displacement, 
but on both the shear slip and normal opening as given in equation 3.3: 
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^'^l (3.3) 

T (u n , U t )) 




Crack tip 

Figure 3.2: Mixed-Mode oscillatory field at crack-tip 



Consequently, the effective traction vector and the corresponding effective displacement can 
respectively be expressed as: 



Cos V Sin V 



(3-4) 



= J(u n ) 2 + u t 2 = = (3.5) 

y x n/ c CosV SinV 



Thus, for a local Mixed-Mode fracture, the critical magnitude of the traction vector now 
depends on the ratio between the shear and normal tractions, which by definition is the phase angle 
given in equation 3.6. 



V = tan 



In effect, as V increases, the normal stress-crack opening curve diminishes from V = 0°, 
while the shear stress-sliding curve expands towards a maximum for W = 90°. 

It should be noted that the phase angle (V) defined here can be at variance from the global 
phase angle (W) defined earlier in equation (2.16b) for the LEIFM, though they both measure the 
relative proportion of the effect of Mode II fracture to Mode I fracture on the interface. For many 
practical systems (Buyukozturk and Hearing, 1998), including cementitious interface, the effect of 
non-zero (/?) is of secondary consequence; hence, the global phase angle can conveniently be 

reduced to V — arctan (~~)as shown earlier. 

Therefore, since the phase angle given by equation 3.6 is of local Mixed-Mode effect, its 
numerical value may vary along the interface - from element to element - (Mei et al, 2010);its 
degree of variation with respect to delamination length is however expected to be relatively 
insignificant for short crack limit where delamination length (Z d ) is less that the overlay thickness 
(^overlay)- Hence, it is appropriate to assume a constant steady-state phase angle during the analysis, 
while treating the interface toughness as independent of the delamination length, but a dependent 
function of the phase angle, so that the interface attains its critical fracture condition when the 
Mixed-Mode energy release rate G ic equals the fracture toughness of the interface G ic (W) as given in 
equation 3.9: 

G ic = G ic QF) (3.9) 
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The expression given here is similar to the one given in equation 2.18, except that the value 
of (f) is local while that of (W) is global. 

4.0 DELAMINATION MODEL FOR DIFFERENTIAL LENGTH CHANGE AT 
DISCONTINUITIES 

Consider a finite bonded concrete overlay system shown in Figure (4.1a) resting on elastic 
foundation, and experiencing a differential length change, either due to thermal gradient or drying 
shrinkage of the overlay. The curling effects of the uniaxial edge-stress condition, say for 
contraction, can be idealized as shown in Figure (4.1b) such that stresses are assumed as acting on a 
beam at the edges rather than on the entire slab surface (Houben, 2006). With increased deformation 
concentrating at the top edge surface of the overlay during the curling process, a partial (l cz ) or true 
interface delamination length (l d ) may be induced along the edges. Apparently, for Mixed-Mode 
delamination to occur, the failure condition given in equation 3.9 must be satisfied. 





blasiichaunaatioi. 



3— f /\ 

b) ^£7nsrlc Fo!ltidGti<^i 



Figure 4.1: Overlay Edge Deformation and Delamination 

From the idealized model shown in Figure 4.1, a 2D plane strain analysis can be implemented 
to estimate the magnitude of the interface delamination driving energy as a function of the overlay 
structural scale and elastic mismatched properties of the bi-material. The structural scale corresponds 
to the ratio of the prescribed overlay thickness to the total thickness of the BCO system in order to 
keep the scale dimensionless. The mismatched elastic properties are controlled by Dundur's 
parameters given in equations 2.10 and 2.11. 

From here, three distinct variables can be associated with the delamination function (D) 
based on the expression given in equation 4.1. 



Hotal 



(4.1) 



From equation 4.1, the first parameter in the bracket ( h ° verla y \ denotesthe structural scale, 

^ h total ' 

and often helps to investigate the thickness response of the overlay to delamination when one of the 
two Dundur's parameters is fixed. In most applications, the limiting value of h overlay falls between 
50 and 125mm. For reasons given earlier, the effect of non-zero (/?) is secondary and subsequently 
can be held fixed for all possible material combinations. Thus, in order to estimate the FPZ (l cz ) 
along the interface, the following relationship holds (Gdoutos, 2005): 



l cz = D GcO = D (^p) 



(4.2) 
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Where, D is as defined in equation 4.1, l ch is Hillerborg's characteristic length defined by 

(E* G \ 
2 Fm ),E* = average elastic modulus of the bi — material (same as equation 2.1S),A m — 

effective traction for mixed mode damage initiation (same as equation 3.4), 
G Fm = Mixed — Mode fracture energy = Gj + G n 



G, = Mode I fracture energy = -X m u^ n cos 2l ¥(A3) 



G n = Mode II fracture energy = -A^u^sin^iAA) 



2 



Um — critical (effective ) mixed mode dis[lacement 

If the Mixed-Mode delamination criterion is specified in terms of fracture energy, equations 
4.5 and 4.6 hold: 



(^L) + (^u) = 1 ( 4. 5) 

\G] C / \Gii c / 

Since, G Fm = G, c + G„ c = \x m u c m = G lc [cos 2 ^ + ^-sin 2 v] (4.6) 

Where, G Fm = fracture toughness or critical fracture enegy for Mixed — Mode 

Gic = Fracture toughness or critical fracture energy in pure mode I 
G llc = Fracture toughness or critical fracture energy in pure mode II 

Rearranging equation 4.2 and expressing the resulting energy release rate in terms of the 
overlay structuralsize, equation 4.7obtains: 

G F m = D (^rMy, a , /?) A " h °7 Hay (4.7) 

With respect to equation 3.9, the delamination failure definition given in equation 4.7 can 
further be expressed as a function of the normalized interface toughness such that: 

r n /-^-overlay 0 \ & Gi c QV) 

G ic = D (- y - ,a,P) = -T- (4.8) 

"■total A m" overlay 

In this respect, the delamination failure coefficient(D) can be numerically estimated as a 
function of the normalized structural scale for different values of (a). This approach relates to plane 
strain problems and similar model has been presented elsewhere (Mei et al, 2010), though with a 
structural scale adjustment. In the literature, several values of (D)based on plane stress problems also 
exist and they are reported in Turon, et. al.(2007). As illustrated in Table 4.1, such values range 
between 0.21 and 1.0; though Hillerborg's and Rice's models where values of (D) approach or equal 
to unityare mostcommon in practice. 

Table 4.1: Cohesive zone length and equivalent delamination dimensionless par ameter 



Model 


^cz 


D 


Hui 


2/3n .E Gc/X cm 2 


0.21 


Irwin 


1/n . E Gc/X cm 2 


0.31 


Dugdale, Barenblatt 


n/8.EGc/X cm 2 


0.40 


Rice, Falk 


9n/32.E Gc/A cm 2 


0.88 


Hillerborg 


E Gcj X cm 


1.00 
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Evidently, from Table 4.1, there exists no unified value of (D)per se. By inspection, the 
values of(D)will depend on the overlay structural scale, the type of problem (plane stress or plane 
strain), the method and the magnitude of loading, and the degree of mismatched elastic properties 
between the overlay and the substrate. 

5.0 CONCLUSIONS 

The research showed that for a composite Bonded Concrete Overlay system, the use of 
numerical computational tools is vital considering the level of complexity involved in determining 
the effects of intrinsic structural and mismatched material properties on interfacial delamination. An 
approach based on plane strain analysis within the context of Interface Cohesive Zone Model has 
been presentedas viable for simulating and predicting delamination mode of failure in BCO systems. 
From the information given in this paper, the following can be concluded: 

• Many of the drawbacks associated with stress-based approach and classical energy-based 
method (LEIFM) such as size effects, incomplete computational analysis and stress 
singularities are overcome in the nonlinear interface fracture mechanics (NLIFM) approach. 

• A unified model where both interfacial crack nucleation and propagation processes are present 
can be implemented using nonlinear interface fracture mechanics approach. 

• The numerical values of delamination failure coefficient(D) and Mixed-Mode energy release 
rates (G ic ) can vary depending on the overlay structural scale, the type of problem (plane stress 
or plane strain), the method and the magnitude of loading, and the degree of elastic 
mismatched properties between the overlay and the substrate. 
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